Molecular Gas and Star Formation in the Nucleus of IC 342: 
C ls O and Millimeter Continuum Imaging 



O 
O 

(N 

C 

a 
in 

(N 



> 



a: 

6 
to 

in 
c3 



David S. Meier, & Jean L. Turner 

Department of Physics and Astronomy, University of California, Los Angeles, CA 90095-1562 
email: meierd@astro.ucla.edu; turner@astro.ucla.edu 



ABSTRACT 

We present high resolution maps (~ 2 ) maps of the J=l-0 and J=2-l transitions of C 18 in 
the central ~ 150 pc of the gas-rich nucleus of IC 342, made with the Owens Valley Millimeter 
Array. From the C 18 maps we are able to obtain the most accurate map of N# 2 to date for IC 
342. Due to their low opacities, the transitions of C 18 give a more reliable estimate of the true 
molecular gas column density distribution than the more common 12 CO and 13 CO isotopomers. 
The morphology of the C 18 emission in the nucleus is a mini-spiral similar to that of the main 
isotopomer, 12 CO, except it is more symmetric, lacking the enhancements to north. We suggest 
that the asymmetries present in 12 CO images may reflect the viewing perspective of the starburst 
region biased by the high optical depths of 12 CO, rather than true asymmetries in the amount 
of molecular gas present. The giant molecular clouds seen in C 18 appear to be nonspherical, 
probably due to tidal arm shearing. Column densities determined from C 18 observations, 1.3 
mm dust continuum, and the virial theorem indicate that the standard Galactic conversion factor, 
Xco, overestimates the amount of molecular gas in the center of IC 342 by a factor of ^2-3 at the 
molecular cloud peaks, and by >3 in the diffuse gas away from the starburst. Revised molecular 
masses based on this conversion factor imply that star formation efficiencies in the starburst 
region are very high. From the distribution of gas and star formation it appears that the sites of 
star formation are dynamically determined, rather than driven by density peaks. Near the central 
star-forming region, evidence is seen for chemical enrichment of C 18 due to massive stars. 
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Our current knowledge about molecular gas in 
galaxies is primarily based on the bright, optically 
thick lines of 12 CO, using the Galactic standard 
conversion factor, Xco- These assumptions, while 
probably adequate when applied to the disks of 
large spiral galaxies such as our own (Solomon et 
al. 1987; Scoville & Sanders 1987; Wilson & 
Scoville 1990; Nakai & Kuno 1995), may not ap- 
ply in the unusual conditions present in nuclear 
starbursts. Enhanced metallicities, densities and 
temperatures, large turbulent linewidths, and in- 
tense radiations fields that are present in galactic 
centers may change the value of Xco ( e g-> Mal- 
oney & Black 1988; Elmegreen 1989; Sakamoto 



1996). It has been suggested that the Xco over- 
estimates gas masses in our own Galactic Center 
by factors of 5 - 16 (Sodroski et al. 1994; Dah- 
men et al. 1998), and in nearby starbursts (eg., 
Smith et al. 1991; Xie, Young & Schloeb 1994; 
Mauersberger et al. 1996; Braine et al. 1997). 

Optically thin lines are promising probes of the 
molecular gas, because they are able to "sample" 
all the column density and not just conditions at 
cloud edges. Unfortunately, because the abun- 
dances of the isotopic species are much lower, the 
emission lines are relatively weak. Only recently 
has instrumentation reached the point where inter- 
ferometeric studies of isotopomers in other galax- 
ies are feasible. High resolution studies of 13 CO 
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(e.g. Hurt & Turner 1991; Turner & Hurt 1992, 
hereafter TH92 ; Papadopoulos, Seaquist & Scov- 
ille 1996; Aalto et al. 1997; Neininger et al. 1998; 
Matsushita et al. 1998; Reynaud & Downes 1999; 
Meier, Turner & Hurt 2000, hereafter MTH00), in- 
dicate that even 13 CO is optically thick in the mas- 
sive clouds of starburst nuclei (Wall & Jaffe 1990; 
Papadopoulos, Seaquist & Scoville 1996; MTH00 
). Therefore, C 18 may be the best hope for deter- 
mining true molecular gas distributions of nuclear 
starburst regions. 

IC 342 provides an excellent galaxy for study- 
ing the relative locations of different molecular gas 
tracers within the nucleus of galaxies because it is 
the closest large spiral with active nuclear star for- 
mation (1.8 Mpc, Table 1: McCall 1989; Madore 
& Frecdman 1992; Karachentsev & Tikhonov 
1993; see Buta & McCall 1999 for a detailed 
discussion of the distance to IC 342). IC 342 is 
nearly face-on, with widespread resolvable molec- 
ular line emission both in dense clouds and a dif- 
fuse medium (eg., Lo et al. 1984; Ishizuki et al. 
1990; Downes et al. 1992; Wright et al. 1993). 

In this paper, the true distribution of the molec- 
ular gas in the nucleus of IC 342 will be investi- 
gated. Questions such as the following will be ad- 
dressed: what is the distribution of H 2 ? What is 
Xco? Does Xco vary over the starburst region? 
and does it correlate with the observed physical 
conditions of the molecular gas? With these goals 
in mind, we obtained ~ 2 resolution Owens Val- 
ley Millimeter Interferometer images of C 18 0(2-l), 
C 18 (l-0), 2.6 mm and 1.3 mm continuum of the 
central region of IC 342, giving spatial resolutions 
of - 15 - 20 pc. The C 18 0(2-l) was observed 
simultaneously with C 18 O(l-0) to allow determi- 
nation of T ex , which is necessary for the determi- 
nation of optically thin column densities. These 
observations are combined with previous data of 
the J=2-l and 1-0 transitions of 12 CO and 13 CO 
(TH92 ; Turner, Hurt & Hudson 1993; MTH00 
) to obtain a detailed and accurate description of 
the molecular gas column density within the cen- 
tral 150 pc. 

2. Observations 

Aperture synthesis observations of both C 18 0(1- 
0) (109.782 GHz) and C 18 0(2-l) (219.560 GHz) 
were made with the Owens Valley Radio Observa- 



tory (OVRO) Millimeter Interferometer between 
1997 October 22 and 1997 December 15 (Table 
2). The interferometer consists of six 10.4 meter 
antennas with cryogenically cooled SIS receivers 
(Padinctal. 1991; Scoville et al. 1994). System 
temperatures (single sideband) ranged from 250 
to 400 K at 2.6 mm and 450 - 1000 K at 1.3 mm. 

The transitions were observed simultaneously 
making use of OVRO's multiple- line spectrome- 
ter capabilities. Therefore, both transitions have 
consistent instrumental configurations, phase cen- 
ters, and weather. The (2-1) data have the same 
(u,v) pattern but with twice the scale as the (1-0) 
data. A nitcrbank with 32 4 MHz channels was 
used, giving velocity resolutions (bandwidths) of 
10.9 km s" 1 (349 km s" 1 ) at 2.6 mm and 5.5 
km s" 1 (175 km s" 1 ) at 1 mm. IC 342's sys- 
temic velocity with respect to LSR is 35 km s _1 
and was centered at channel 16.5. Two pointings 
with phase centers of (51950) = 03' l 41 m 57. s 9; 
Si (£1950) = +67°56'29."0 and a 2 (B1950) = 
03 h 41 m 57. s 0; <5 2 (B1950) = +67°56'26."0 were 
observed to cover the entire nuclear mini-spiral. 
Data were calibrated using the MMA software 
package. Calibrators 3C84 and 0224+671 were 
observed in 20 minute intervals for phase calibra- 
tion. Absolute flux calibration was done using 
Neptune and Uranus as primary flux calibrators 
and 3C273 as a secondary flux calibrator for both 
transitions. Absolute flux calibrations are esti- 
mated to be good to 10 - 15% for the 2.6 mm 
data and 20 - 25% for the 1.3 mm data. 

The dataset was mosaicked using the MIRIAD 
package. The maps are naturally weighted and 
primary-beam corrected. Since the maps are mo- 
saics, noise levels vary across the maps. Noise 
levels reported throughout this paper are those 
measured from line free channels of the spectral 
cube half-way between the map center and the 
FWHM points. The true noise level is slightly 
lower (~ 10%) at the phase centers and somewhat 
higher towards the edge of the primary beams. All 
further data reduction was done with the NRAO 
AIPS package. In making the zeroth moment (in- 
tegrated intensity) maps, only emission greater 
than 1.2cr was included. Due to the fact that 
the signal-to-noisc (S/N) is lower in the C 18 0(2-l) 
map, the choice of 1.2a could, in principle, lower 
the integrated intensities of the C 18 0(2-l) map 
relative to the higher sensitivity C 18 O(l-0) map. 
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However, the ratio of the integrated intensity maps 
made with no clipping is identical to those with 
clipping to within 5%, indicating this effect is not 
important. 

The (u,v) coverage for these observations imply 
that emission on scales larger than —36 (3 mm) 
or -18 (1.3 mm) is resolved out. To estimate 
the amount of extended flux the interferometer re- 
solves out in the C 18 O(l-0) map, it was compared 
with the single-dish spectra from the 30 m IRAM 
telescope (Eckart ct al. 1990) at 21 resolution. 
The integrated intensity of the OVRO map when 
convolved to 21 is 3.6 K km s _1 , which is —67% of 
the 4.8 K km s -1 detected by IRAM. This is con- 
sistent with the 70 - 75% detected in the 12 CO(1-0) 
map (Levine et al. 1994b), and not much lower 
than the 80% estimated for the 13 CO transitions 
(TH92 ; MTH00 ). A similar analysis yields 4.9 
K km s- 1 for C 18 0(2-l) or -60% of the single- 
dish flux. Therefore, the amount of flux resolved 
out in 12 CO(1-0), 13 CO, C 18 are all similar and 
should not strongly effect the inter-comparisons. 

3. Results 

3.1. The C ls O Morphology 

The C 18 O(l-0) and C 18 0(2-l) channel maps 
are presented in Figures 1-2. Peak antenna tem- 
peratures obtained in the 4 MHz channels are — 
0.75 K for C 18 O(l-0) and 1.5 K for C 18 0(2-l). 
Emission is present in channels from Vlsr = -14 
km s _1 to 73 km s _1 with higher velocities in the 
northeast. The line center and velocity dispersions 
of both transitions of C 18 are consistent with 
those of 13 CO (TH92 ). Integrated intensity maps 
of the C 18 O(l-0) and C 18 0(2-l) transitions, along 
with the 12 CO(1-0) map (Levine ct al. 1994b) are 
plotted to the same scale in Figure 3. Integrated 
intensities of each transition are listed in Table 3 
for various locations across the nucleus (nomencla- 
ture from Downes et al. 1992, Figure 2 of MTH00 

)■ 

The basic morphology of the nuclear mini-spiral 
seen in C 18 O(l-0) is similar to the 12 CO(1-0) tran- 
sition, but there are noticeable differences. Ob- 
servations of 12 CO and HCN show that there are 
five major peaks of molecular gas in the nucleus 
(Ishizuki et al. 1990; Downes et al. 1992). GMC 
C, along the northern arm, is brightest in 12 CO, 
GMC B is the molecular cloud associated with 



the starburst, and GMC A is the nearest to the 
dynamical center of IC 342. GMCs D and E are 
more distant GMCs along the northern and south- 
ern arms, respectively. All GMCs except GMC A 
are readily seen in C 18 O(l-0). GMC A is substan- 
tially weaker in the 1-0 transition of C 18 relative 
to the 12 CO(1-0) transition. However, in the 2-1 
transition of C 18 0, GMC A is similar in relative 
strength. The weakness of GMC A in C 18 O(l-0) 
is probably a temperature effect (§3.2.1). 

Another difference is that C 18 O(l-0) distribu- 
tion, like 13 CO(1-0) (TH92 ; Wright et al. 1993), 
appears more symmetric than 12 CO (1-0). In 
12 CO(1-0), the northern arm is brighter than the 
southern arm by about 50% (the integrated flux 
of the entire northern arm is 1200 Jy km s _1 , 
while the entire southern arm is 800 Jy km s _1 ). 
In C 18 O(l-0) it is only -15% percent brighter 
(the integrated flux of the entire northern arm 
is 15.4±2.5 Jy km s _1 , as compared to 13.6±2 
Jy km s _1 for the southern arms). C 18 O(l-0) 
preferentially arises from the outside (leading) 
edges of the arms, similar to the 13 CO(1-0), while 
12 CO(1-0) comes from broad arms with abundant 
diffuse emission upstream from the arms. The 
less abundant isotopomers do not show the broad 
diffuse arms nearly as dramatically as the 12 CO 
(Wright et al. 1993). Although the sensitivity of 
the C 18 O(l-0) maps are not extremely high, the 
high sensitivity 13 CO(1-0) line displays the same 
trend. These differences between 12 CO and the 
rarer isotopes can be understood if the geome- 
try and excitation of the nuclear gas clouds are 
considered. 

3.1.1. The Geometry of the Starburst 

From Figure 4, it can be seen that the two 
molecular arms in the nucleus of IC 342 progress 
inward and wrap around the central starburst. 
The arms appear to terminate in a partial ring 
surrounding the burst. If the two arms are due to 
a trailing density waves with orbital motions pre- 
dicted for barred spiral arms outside of the Inner 
Lindblad Resonance (ILR) (eg., Roberts, Huntley 
& van Albada 1979; Athanassoula 1992), then the 
velocity pattern implies the northern arm is be- 
hind the starburst in projection and the southern 
arm is in front. Comparisons of Ha to thermal ra- 
dio continuum and observations of B7/Ha confirm 
that extinction towards the southern portion of 
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the starburst complex is substantially higher than 
towards the north (Figure 4; Figure 5 of MTHOO 
; Boker et al. 1997). 

This geometry implies that from our perspec- 
tive, we are seeing the "starburst lit" face of the 
molecular clouds when viewing the northern arm 
(seeing the surface PDR regions), and the "back- 
side" of the molecular clouds for the southern 
arm (seeing the ambient molecular gas). For op- 
tically thick transitions such as those of 12 CO, 
which trace the cloud photospheres (Turner, Hurt 
& Hudson 1993), the lines should appear brighter 
along the northern arm relative to the southern 
arm, hence are sensitive to temperature effects. 
Whereas the optically thin transitions will trace 
column density. Therefore this asymmetric pat- 
tern can manifest itself in optically thick transi- 
tions, but not in transitions that are either inher- 
ently optically thin due to opacity, or effectively 
thin due to large velocity gradients. 

3.1.2. Giant Molecular Clouds in IC 342 

Because of its low opacity and higher critical 
density, C 18 0(2-l) emission, like 13 CO(2-l), is 
dominated by the dense giant molecular clouds 
(GMCs) seen in HCN emission (Downes et al. 
1992). The C 18 0(2-l) observations have the high- 
est spatial resolution yet obtained (^13 pc) in 
IC 342, sufficient to resolve the nuclear GMCs 
into substructure. Therefore, C 18 0(2-l) is used 
to derive the properties of the individual GMCs. 
Each cloud, identified only as a region of spatially 
and spectrally localized emission, not necessarily a 
gravitationally bound entity, has been fit with an 
elliptical Gaussian to determine its size, linewidth, 
location and intensity. Of the five GMCs identi- 
fied by Downes et al. , only GMCs A, B, C and E 
are fit since GMC D is at the edge of the C 18 0(2- 
1) primary beam. The four major GMCs are re- 
solved into eight components (numbers are added 
to the letter designations in the case of substruc- 
ture) . The results of the fits are tabulated in Table 
4. 

All clouds are resolved (fitted sizes >l/2 the 
FWHM of the beamsize) except cloud C3. The 
GMCs are composed of several smaller clouds that 
are ~ 10 x 30 pc in size, with masses of about 
~ 10 6 M and FWHM linewidths of 20 - 30 
K km s _1 . The GMCs tend to be non-spherical 
with axial ratios ~2 and position angles gener- 



ally oriented in the direction of the spiral arms. 
This is consistent with the presence of arm stream- 
ing motions over the nuclear region (Downes et 
al. 1992; TH92 ). Beam averaged densities for 
the molecular clouds that make up the GMCs are 
<~ 10 3 ' 5 cm~ 2 (§3.2.3). These values are quite sim- 
ilar to those found for the "envelope" of the Galac- 
tic Center GMC, Sgr B2 (eg., Gordon et al. 1993; 
Hiittcmcister et al. 1995). 

3.2. Comparing CO Isotopomers: Gas Ex- 
citation and Opacity in IC 342 

3.2.1. Excitation Temperatures from C 18 0(2- 
l)/C ls O(l-0) 

Comparisons between rotational lines with dif- 
ferent J states of optically thin transitions allow 
for the mapping of gas excitation temperatures 
(T e x)- Generally, the main isotopomers of CO 
are optically thick in galactic nuclei and as such 
are of little use in constraining T ex (eg., Braine 
& Combes 1992). Moreover, gradients in tem- 
perature over small scales appear to dominate the 
12 CO line ratios rather than T ex (Turner, Hurt 
& Hudson 1993). Optically thin line ratios, on 
the other hand, can be quite constraining, mak- 
ing the C 18 O(2-l)/C 18 O(l-0) line ratio a sensitive 
measure of T ex in IC 342. The ratio of the C 18 0(2- 
1) integrated intensity to the C 18 O(l-0) intensity, 
under the LTE assumption is: 

/ 8 T 2 idv ^ 
J 18 T w dv ~ 

18 f2i( 18 J2i(T ex ) - 18 J 21 (T cmfc ))(l - e- 18 ^) 
18 /io( 18 Jio(T e:c ) - 18 Jio(T cmfc ))(l - e- 18 -)' 

where J u (T ex ) = (hu / k) / (exp{hu / kT ex } - 1), r 
and / are the optical depth and filling factor of 
each transition, respectively. 

In Figure 5a, the C 18 O(2-l)/C 18 O(l-0) line ra- 
tio map is presented for the nucleus of IC 342. 
Ratios of integrated intensities are used instead 
of the peak main-beam temperatures to increase 
the S/N. This requires that the C 18 0(2-l) and 
C 18 O(l-0) line profiles be the same. This is a 
good approximation (compare Figure 1 and 2). In 
generating the ratio maps, the integrated intensity 
maps for C 18 0(2-l) and C 18 O(l-0) were convolved 
to matching beamsizes, divided, and blanked any- 
where where emission was < 3tr in either map. 
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Errors in the ratio maps are conservatively esti- 
mated to be <35% in magnitude and < 1 in po- 
sition. 

Aside from the two peaks seen in the ratio map, 
the C 18 O(2-l)/C 18 O(l-0) ratio map has relatively 
constant values of 0.7 - 1.3 over much of the cen- 
tral mapped region, with lower values preferen- 
tially found in the south. The highest C 18 0(2- 
l)/C 18 O(l-0) peak, with a value of 2.2, is found 
at the starburst region. The heating source for 
the warm molecular gas at this location is proba- 
bly the massive stars of the starburst. The second 
ratio peak is at GMC A. GMC A is unusual in that 
it is not near either of the major radio continuum 
sources (§3.3), which suggests that the heating of 
GMC A is not due to star formation. Some dy- 
namical heating mechanism may be at work here, 
such as shock heating due to large scale bar-driven 
gas flows, which is consistent with the fact that it 
is the GMC closest to the dynamical center of the 
galaxy, and to the ring of hot shocked gas seen in 
H 2 (Bokcr et al. 1997). 

A map of T ex for the nuclear region derived us- 
ing the C 18 (2-l)/C 18 (l-0) ratio is displayed in Fig- 
ure 5b. T ex ranges from 7 - 10 K over much of the 
mapped region. The T ex of the gas near the star- 
burst region is 19±8 K. The secondary T ex peak 
at GMC A has T ex = 13±4 K. Excitation temper- 
atures derived from the C 18 transitions are lower 
than those estimated from 13 CO (MTH00 ). We 
estimate filling factors using 18 /io ~ T m {,/T e:E Tio. 
For T ex ~ 10 K and t w ~ 0.1 - 0.2 (§3.2.4), 
the implied filling factors are 18 /io ~ 0.5 — 1, 
which given the different beamsizes, are consistent 
with those found for 13 CO. This implies that the 
C 18 0(2-1) emission detected arises predominately 
from the same molecular gas component as the 
13 CO(2-l) emission (MTHOO ). 

3.2.2. Molecular Gas Opacities and Abundances: 

12 CO (1-0) /G 18 O(l-0) & 13 GO (1-0) /G 18 0(1- 
0) 

Comparison of the abundant and optically thick 
12 CO line with its less abundant and thinner 
cousins, 13 CO and C 18 0, can be used to trace 
optical depths and place limits on isotopic abun- 
dances. The isotopic line ratios are displayed in 
Figure 6a-c and Table 3. All of the isotopic line 
ratios discussed in this paper are between the (1-0) 
transitions of each isotopomer and so we suppress 



the (1-0) notation for the remainder of the pa- 
per. The 12 CO/C 18 line ratio ranges from 11 - 
>110 over the mapped region. Regions that are 
blanked due to low S/N tend to correspond to 
high 12 CO/C 18 ratios since the weakness of the 
C 18 O(l-0) line is the limiting factor. At the sites 
of the GMCs, 12 CO/C 18 ratios of 30 - 60 are 
found, while values for 13 CO/C 18 are 5-6. Typi- 
cal Galactic Center values for the 12 CO/C 18 and 
13 CO/C 18 are - 60 - 80 and - 10, respectively 
(eg., Oka et al. 1998; Dahmen et al. 1998). 
As in the Galactic Center (Dahmen et al. 1998), 
variations in this ratio in the nucleus of IC 342 are 
substantial over small (~30 pc) spatial scales. The 
highest column density regions at peaks B, C, and 
E tend to have the lowest isotopic ratios. This is 
expected since these are presumably the optically 
thickest regions. 

The Central Trough region has surprisingly low 
values of 12 CO/C 18 and 13 CO/C 18 line ratios. 
The Central Trough region, which is a site of weak 
emission in all observed transitions, is the only 
firmly established region of optically thin 13 CO gas 
over the central distribution (MTHOO ). As such, 
it is expected to have high isotopic ratios and give 
the actual 13 CO/C 18 isotopic abundance ratios. 
Surprisingly, both isotopic ratios, 12 CO/C 18 = 
24 and 13 CO/C 18 = 3.2, are lower here than 
the average across the rest of the spiral arms 
(C 18 over luminous). Along the spiral arms, the 
12 CO tends to dominate C 18 (highest ratios) on 
the upstream sides of the arms. C 18 becomes 
relatively more prominent as one moves towards 
the leading edge of the spiral arms (low ratios). 
The same effect is seen in 12 CO/ 13 CO observa- 
tions with higher sensitivity to extended structure 
(Wright et al. 1993). 

Much of the structure seen in the isotopic ra- 
tio maps appear to be due to excitation effects. 
The peak in the 12 CO/C 18 ratio map at GMC 
A and towards the starburst near GMC B result 
from depopulation of the 1-0 level in the opti- 
cally thin C 18 due to the increased T e x, as con- 
firmed by the high C 18 O(2-l)/C 18 O(l-0) ratios. 
The trend of decreasing ratio counter-clockwise 
across the arms likely results from changes in den- 
sity across the spiral arm due to the density wave 
(Wright et al. 1993). The variation in excitation 
demonstrates the importance of observing higher J 
transitions when dealing with optically thin tran- 
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sitions. If one relies, as is often the case for ex- 
tragalactic studies, on just the observations of the 
(1-0) transition of CO and its isotopomers, it is 
possible that variations in excitation can be mis- 
taken for variations in morphological structure or 
column density. 

A lower limit to the true 12 CO/C 18 abun- 
dance ratio can be estimated from the observed 
ratio. In the limit of both transitions being opti- 
cally thin (LTE), the isotopic line ratio will re- 
flect the true isotopic abundance ratio. Since 
12 CO(1-0) is not optically thin, the ratio is a lower 
limit. For conservative lower limits to [ 12 C/ 13 C] 
and [ 16 0/ 18 0], we take the highest (trustworthy) 
values of the ratio. The highest observed value of 
the 12 CO/C 18 line ratio is -110. The [ 16 0/ 18 0] 
abundance ratio must be higher than this. The 
peak values of the 13 CO/C 18 line ratio are about 
11. However, both of these values are overesti- 
mates because of the non-LTE effects discussed 
in §3.2.4. When we have corrected for the non- 
LTE effects, we estimate [ 16 0/ 18 0] > 110 and 
[ 12 C/ 13 C] > 20. Though not well constrained, 
the values are consistent with those assumed as 
well as the (also weakly constrained) values previ- 
ously estimated for IC 342 based on other molec- 
ular species (Hcnkcl ct al. 1998). 

The low isotopic ratios seen towards the Cen- 
tral Trough appear to require truly anomalous 
abundances. The low ratios in this region imply 
that C 18 is overabundant relative to 12 CO and 
13 CO. Neither isotope-selective photodissociation 
nor isotope-selective fractionation can be the ex- 
planation because in both cases the abundance of 
12 CO and 13 CO are predicted to be larger relative 
to C 18 0, not smaller (Langcr et al. 1980, 1984; 
van Dishoeck & Black 1988). This makes the most 
likely possibility for the higher abundance of C 18 
enrichment from 18 O rich, massive star ejecta. A 
complete knowledge of the cooking sites of 18 O has 
yet to be achieved, but it appears that 18 O is a 
secondary product of massive stars (see eg., Wil- 
son & Matteucci 1992; Henkel & Mauersberger 
1993; Wilson & Rood 1994; Prantzos, Aubert & 
Audouze 1996). If the central starburst is bi- 
ased towards massive stars, a relative increase in 
the abundance of 18 (and hence C 18 0) is pos- 
sible. Low 12 CO/C 18 values have been found 
before in other starburst galaxies such as NGC 
253, M 82, NGC 4945 as well as the AGN NGC 



1068 (eg., Henkel & Mauersberger 1993; Mauers- 
berger et al. 1996; Papadopoulos, Seaquist & 
Scoville 1996). Using a chemical evolution model, 
Henkel & Mauersberger (1993) have modeled the 
isotopic abundance ratio for a starburst typical of 
a galactic nuclei like M 82, NGC 253 and possibly 
IC 342. They find that the [ 16 0/ 18 0] can drop 
as low as —150. IC 342 also has a rather high 
C 18 O(l-0)/C 17 O(l-0) line ratio of 6.5, consistent 
with enhanced 18 (L. Sage, from Henkel et al. 
1994). Therefore, for the Central Trough region, 
it is possible that the [C ls O/H 2 ] abundance may 
be higher than has been assumed. 

3.2.3. Molecular Gas Density: C 18 O LVG mod- 
eling 

A sample of Large Velocity Gradient (LVG) ra- 
diative transfer models were run to relate the ob- 
served intensities and line ratios to the physical 
properties of the cloud, nn 2 and T& (eg., Goldrc- 
ich & Kwan 1974; Scoville & Solomon 1974; dc 
Jong ct al. 1975). LVG models are probably not 
appropriate for all lines in these tidally disturbed 
clouds, (Turner, Hurt & Hudson 1993; MTH00 ), 
but arc nonetheless instructive, particularly when 
the transitions get optically thinner, as is the case 
for 13 CO and C 18 0. The LVG model used is the 
same as discussed in MTH00 , except a [C 18 0/H 2 ] 
abundance of 3.4 x 10~ 7 is used (§4.1). The ratio 
of cloud lincwidth to cloud size is used to con- 
strain the velocity gradient for the models (Table 
4). Typical values of the velocity gradients for the 
GMCs are ~ 1 — 5 km s _1 pc _1 . In Figure 7, we 
display two models with velocity gradients of — 1 
and 5 km s _1 pc _1 , for both 13 CO and C 18 (see 
MTH00 , for details of the 13 CO models). 

Since C 18 is optically thinner than 13 CO, it is 
less affected by radiative trapping and thermalizes 
at a slightly higher density than the correspond- 
ing lines of 13 CO. LVG modeling of C 18 may 
be capable of constraining the densities, whereas 
13 CO has only provided lower limits to the den- 
sity. Best fit solutions obtained from the C 18 
temperatures and line ratios are indeed close to 
the lower end of the solution range found in 13 CO, 
and well constrained. The C 18 observations 
imply that most of the central GMCs have < 
n H2 > ~ 10 3 " 3 - 5 cm" 3 and T fc ~ 10 - 40 
K. These moderate densities come directly from 
the fact that T ex derived from C 18 are lower 
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than those derived from 13 CO. These densities are 
slightly smaller than those implied from the HCN 
observations (Downes ct al. 1992), probably in- 
dicating that there is some further clumping that 
is not resolved. Closer inspection shows that on 
average the C 18 solutions for GMC A, B and 
C imply slightly warmer, lower density gas than 
does 13 CO. This is probably caused by the fact 
that 13 CO and C 18 were assumed to have the 
same filling factor (0.25 over a 4. 8x3. 9 beam; 
MTH00 ). Since the critical density is higher for 
C 18 0, one might expect cloud sizes to be slightly 
smaller for C 18 than for 13 CO. If the relative fill- 
ing factor for the C 18 emitting gas is lowered, the 
models become consistent for all locations. In ev- 
ery case, the decrease in filling factor required for 
agreement between 13 CO and C 18 are less than a 
factor of two. It is interesting to compare the aver- 
age densities derived from the LVG analysis with 
what is predicted from high resolution C 18 0(2-l) 
and the virial theorem (§4.1). The densities esti- 
mated from both methods are not too different, 
with the values derived from the virial theorem 
tending to be about 1.5-3 times larger than the 
LVG models. 

3.2-4- Isotopic Line Ratios: Complications and 
Important Considerations 

These multi-line, multi-transition data on IC 
342 provide a wealth of information on the com- 
plex interactions between opacity and temperature 
for molecular clouds in starburst environments. 
Here we consider what these multi-transition stud- 
ies as a whole indicate for the opacities of the CO 
lines. The magnitude and effect of several differ- 
ent non-LTE effects are estimated for the nucleus 
of IC 342. 

For IC 342, the T ex and filling factors derived 
from the 12 CO, 13 CO and C 18 lines can depend 
on the transition considered. A single T ex or fill- 
ing factor cannot be assumed for all the lines. 
This complicates the determination of opacities 
and abundances from the line ratios. Since 13 CO 
and C 18 have lower opacities, their effective crit- 
ical densities are higher than for 12 CO. They are 
more confined to the cloud cores and will have 
smaller relative filling factors (Wall & Jaffe 1990; 
MTH00 ). For the externally heated clouds in 
IC 342, the 12 CO transitions will also have higher 
T ex than the optically thin species (Turner, Hurt 



& Hudson 1993; MTH00 ). Therefore the iso- 
topic ratios between 12 CO and other isotopomers 
assuming LTE will underestimate the line opaci- 
ties. 

It is difficult to estimate the degree to which 
differing filling factors impact the line ratios in- 
dependent of the opacity. However, differences in 
the T ex can be estimated since T ex have been sep- 
arately derived for 12 CO, 13 CO and C 18 (Turner, 
Hurt & Hudson 1993; MTH00 , §3.2.1). If 
the peak brightness temperature of the 12 CO(2-l) 
transition is taken as a lower limit to the T ea; , 
the 13 CO(2-1)/ 13 CO(1-0) line ratio under the op- 
tically thin assumption to estimate the 13 T ex , 
and similarily the C 18 O(2-l)/C 18 O(l-0) line ra- 
tio for the 18 T ex , then the magnitude of the tem- 
perature differences can be estimated. Over the 
central GMCs, typical derived temperatures are 
12 T ex > 26 K, 13 T ex ~ 20 K and 18 T ex ~ 10 K. 
This implies that the line ratios involving C 18 
are ^2-3 times larger than expected assum- 
ing the same excitation temperature applies for 
all transitions. Taking into account these correc- 
tions, the isotopic line ratios imply 13 ti <~ 0.7 — 3 
( 18 rio ~ 0.1 — 0.5). These results are consis- 
tent with 13 ti > 1 as was previously found from 
the 13 CO line ratios (MTH00 ). When one com- 
pares isotopic line ratios in IC 342 to find abun- 
dances and opacities, the derived values are un- 
derestimates because the optically thinner tran- 
sitions have lower T ex than the thicker transi- 
tions. We expect that this variation in T ex for 
CO lines of different opacities may be common in 
other starburst nuclei, and should be taken into 
account when making estimates of opacities and 
abundances. 

3.3. Millimeter Continuum Morphology 

Continuum emission was detected at 2.6 mm 
and 1.3 mm in IC 342, in a 1 GHz continuum chan- 
nel. The 2.6 mm and 1.3 mm continuum maps, 
convolved to the same beamsize, are displayed in 
Figures 8a-b. The 2.6 mm continuum emission 
shows a twin-peaked structure similar in morphol- 
ogy to the 2 cm and 3.4 mm continuum emission 
(Turner & Ho 1983; Downes et al. 1992). The two 
strongest sites of millimeter continuum are coinci- 
dent with the nuclear starburst (near GMC B) and 
a secondary cluster centered on GMC C, respec- 
tively. The peak intensity of 2.6 mm continuum is 
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4.4 mjy bm . 

Figures 8c-d display the spectral index between 
2 cm/ 2. 6 mm and 2.6 mm/1 mm over the central 
region. The 2 cm to 2.6 mm spectral index ranges 
from a% 6 ~ —0.4 to 0.3 (S v cx v a ) over the map, 
similar to what is found for the spectral index be- 
tween 6 cm and 2 cm (Turner & Ho 1983). The 
generally flat spectrum over the nuclear region is 
consistent with thermal bremsstrahlung emission. 
Lower values of a 2 (~ —0.4) towards the north- 
west side of the central ring imply some contam- 
ination by non-thermal synchrotron emission at 
2 cm in this region. Towards GMC C, a| 6 is 
^+0.2 to +0.3. Assuming all the 2.6 mm contin- 
uum emission is from bremsstrahlung, the implied 
number of ionizing photons, N uv , present at each 
peak is ~ 1.8 x 10 51 s _1 , or about 80 06 stars 
(assuming an T e = 8000 K; eg. Mczger & Hender- 
son 1967; Vacca, Garmany & Shull 1996). The 
total flux in the mapped region is S2.6mm = 26 
mjy, implying that the entire central star-forming 
region (~ 75 pc) has N uv ~ 10 52 s" 1 (~ 460 06 
stars; consistent with that estimated by Becklin 
et al. 1980; Turner & Ho 1983; Downes ct al. 
1992). 

1.3 mm continuum emission is detected only 
from the molecular column density peaks as traced 
by C 18 because the S/N is lower at this fre- 
quency. Peak 1.3 mm intensities arc 8 mjy bm -1 
(~ 4cr) and 12 mjy bm" 1 (~ 6a) for GMC B and 
C, respectively. Both 1.3 mm continuum sources 
have rising spectral indices of a\§ ~+2.5 between 
2.6 mm and 1.3 mm. Rising spectral indices at 
1.3 mm are seen tentatively towards GMC A as 
well (Figure 8d). The 1.3 mm continuum appears 
to be largely thermal dust emission. In fact for 
GMC C, the slightly rising spectral index between 
2 cm and 2.6 mm indicates that a small amount 
of dust emission may even contribute to 2.6 mm 
continuum at this location (< 1 mjy bm -1 ). We 
discuss the dust masses implied by the 1.3 mm 
continuum in the following section. 

4. Discussion 

In order to study the relation of galactic dynam- 
ics, molecular gas, and star formation, one needs 
accurate surface density maps of molecular gas. 
Estimates of column densities made using 12 CO 
and a conversion factor, Xco ar e commonplace, 



but require Xco to be valid, which is not yet es- 
tablished. Comparisons between masses estimated 
using various methods need to be made in order 
to investigate the validity of Xco in starburst en- 
vironments. Given the wealth of data available for 
IC 342, and in particular the valuable C 18 and 
1.3 mm maps, accurate values of molecular column 
density can be obtained. 

4.1. H2 Gas Column Densities: Compar- 
ing C ls O, Xco, Virial Masses and 
Dust 

Since C 18 is optically thin, H 2 column den- 
sities can be obtained from the total C 18 O(l-0) 
intensity: 

N(H 2 )ci SO = 2.42 x 10 14 cm" 2 

5-27 

x — Ia&o {K km s ). 

{eTe* -1) 

We adopt an abundance ratio of [H 2 ]/[C 18 0]=2.94x 
10 6 or [C 16 O]/[C 18 O]=250 (Frcrking, Langer & 
Wilson 1982; Hcnkel et al. 1994; Wilson & Rood 
1994). Figures 4 and 5c shows a map of Nh 2 
calculated from optically thin C 18 0. Within the 
2-1 primary beam, the derived T ex obtained from 
the C 18 O(2-l)/C 18 O(l-0) line ratios are used in 
calculating the column density. Since T ex was 
determined to be - 7 - 10 K (§3.2.1), T ex = 10 
K is assumed when deriving column densities for 
regions outside the C 18 0(2-l) primary beam. Col- 
umn densities measured in the central region of 
IC 342 range from < 0.25 - 4.8 x 10 22 cm" 2 (cor- 
rected for resolved out flux assuming a uniform 
distribution). These values obtained for optically 
thin C 18 can be compared with column densi- 
ties estimated using the a Galactic disk Xco of 
~ 2 x 10 20 cm" 2 (K km s" 1 )" 1 (Strong et al. 
1988; Hunter et al. 1997). Using A co , we ob- 
tain Njj 2 column densities ranging from 2-12 
x 10 22 cm -2 over the ~ 25 pc beamsizes, a factor 
of 3 to 13 higher than the column densities based 
on C 18 0. 

Two additional column density estimates are 
considered. First, since the GMCs in the center 
of IC 342 are resolved, the total molecular mass 
can be estimated from the virial theorem. For 
molecular clouds with density profiles given by 
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p oc R 7 , the virial mass is given by (eg., Ma- 
cLaren, Richardson & Wolfendale 1988): 

^ v ( 3 — 7 ) (fcm s _1 ) (pc) 

where R is the radius of the cloud, and Avi/ 2 is the 
FWHM of the linewidth. The radius of the clouds 
is assumed to be 1.4i?i/ 2 , or 0.7\/~ab (eg., Wilson 
& Scoville 1990; Goldsmith, Bergin & Lis 1997), 
where a and b are the FWHM fitted sizes of the 
major and minor axes. A 7 = 1 is assumed for the 
GMCs, consistent with what is found for Galactic 
clouds (eg., Cernicharo, Bachillcr & Duvert 1985; 
Scoville et al. 1987; MacLarcn, Richardson & 
Wolfendale 1988). However, the derived masses 
are only weakly dependent on the choice of 7. The 
masses will change from their 7 = 1 values by less 
than a factor of 50% for any 7 between 0-2. Virial 
cloud masses for 7 = 1 are displayed in Table 4. 

Gas masses may also be estimated from opti- 
cally thin thermal dust continuum emission. The 
2.6 mm and 1.3 mm continuum maps are used to 
measure the amount of dust. Assuming the clouds 
are entirely molecular and have a constant gas to 
dust ratio of 100 by mass, the gas mass is related 
to the 1.3 mm dust continuum flux by (eg., Hildc- 
brand 1983): 




where k v is the dust absorption coefficient at 
this frequency, S\ mm is the 1.3 mm flux, D is 
the distance and T d is the dust temperature. 
The dust opacity, n v , at 1.3 mm is taken to be 
3.1 x 10~ 3 cm 2 g _1 , but is uncertain by at least a 
factor of four (Pollack et al. 1994). A T d = 42 K 
is assumed based on the FIR colors (Bccklin et al. 
1980; Rickard & Harvey 1984). The dust temper- 
ature applicable to the 1.3 mm observations could 
be lower than this if a cool dust component missed 
by the FIR observations exists. Since dust emis- 
sion is only detected at the two sites with major 
star formation, we do not expect that a substantial 
cool dust component exists. Therefore, uncertain- 
ties in the derived gas mass are dominated by the 
uncertainty in opacity. 



The amount of dust emission is estimated by 
assuming the 2.6 mm continuum emission is en- 
tirely thermal free-free emission. The nearly iden- 
tical morphology and flat spectral index between 
2.6 mm and 2 cm makes this a reasonable assump- 
tion (Turner & Ho 1983). (In the youngest star- 
forming regions, there are regions of compact op- 
tically thick HII regions that do have rising spec- 
tral indices, [eg., Beck et al. 1996; Turner, Ho 
& Beck 1998; Beck, Turner & Kovo 2000], but 
even these regions will be optically thin by 3 mm.) 
Thermal free- free emission at 1.3 mm is removed 
by extrapolating the 2.6 mm map to 1.3 mm us- 
ing a = — 0.1. The remaining flux of ^8 mJy 
bm- 1 at GMC C and -4 mJy bm" 1 at GMC B 
is estimated to be dust emission. Beam averaged 
column densities based on dust emission are shown 
in Table 5. The gas masses implied by the dust are 
7.4 x 10 5 M and 3.6 x 10 5 M (or 7.7 x 10 22 cm" 2 
and 3.8 x 10 22 cm' 2 ), for the two clouds detected 
at 1.3 mm, respectively. Uncertainties are at least 
a factor of 2 - 4, dominated by uncertainties in k v . 
Given the large uncertainties, dust masses should 
be treated as indicative only. Dust emission in 
IC 342 will be discussed further in a later paper 
(Meier et al., in preparation). 

4.2. The Validity of X CO in IC 342's Nu- 
clear GMCs 

The different methods of obtaining M# 2 and 
Nff 2 can now be compared to determine the valid- 
ity of the Xco in the nucleus of IC 342. Column 
densities derived for the GMCs in IC 342 from 
the three different methods described above are 
listed in Table 5. At all locations, the lowest col- 
umn densities are estimated from the C 18 data. 
The C 18 column densities at the GMCs are 3 - 
4 times lower than those found using the Galactic 
Xco- Column densities based on dust emission 
are slightly higher but consistent with those esti- 
mated from C 18 0. Column densities based on the 
virial theorem tend to be the largest of the three 
methods and are approximately equal to the val- 
ues obtained using the Galactic Xco- 

The dominant uncertainties in the column den- 
sity estimates based on the C 18 data arc the 
[C 18 0/H2] abundance ratio and differing source 
sizes. A Galactic Center value of [C 18 0/H 2 ] has 
been assumed. The isotopic line ratios, while not 
strongly constrained, indicate that this value is 



9 



reasonable (§3.2.2). For the C 18 column density 
to be consistent with Xco value, an abundance 
of [ 16 0/ 18 0] ~ 1000 would be required. This 
value is larger than the highest ISM values seen 
in the Galaxy (—600) and at least a factor of four 
larger than is typically seen in galactic nuclei (eg., 
Hcnkel & Mauersberger 1993; Wilson & Rood 
1994). If anything, the nucleus of IC 342 could 
have a [C 18 0/H 2 ] abundance that is higher than 
the 3.4 x 10~ 7 assumed, since nuclear processing 
in massive stars tend to raise the [C ls O/H 2 ] abun- 
dance rather than lower it (§3.2.2). From an abun- 
dance standpoint, we expect that the C 18 values 
listed in Table 5 do not underestimate the molecu- 
lar gas by a large amount. From a photochemical 
standpoint, it is possible, given isotope selective 
photodissociation and gas density variations that 
C 18 emission will remain undetected in regions 
with detectable 12 CO. In IC 342, the cloud sizes do 
appear to be smaller in C 18 than in 12 CO (§3.2.3; 
MTH00 ), so the column densities estimated by 
C 18 cloud be slight underestimates. Gas excita- 
tion can also influence the optically thin estimates, 
although our dual line observations mostly correct 
for excitation. 

Column densities estimated from the virial the- 
orem are very likely overestimates of the true 
molecular gas column density. The clouds tend 
to be somewhat filamentary in appearance, con- 
sistent with tidal shear or arm streaming motions 
(TH92 ; Downcs ct al. 1992). The lincwidths of 
the GMCs are probably larger than virial and the 
clouds may not be bound. In addition, lincwidths 
in this region are slightly overestimated simply by 
the finite resolution of the telescope in the steep 
part of the rotation curve. Not surprisingly, there- 
fore, the LVG models predict average densities a 
factor of 1.5 - 3 times lower than the virial predic- 
tion. 

While the dust masses are more uncertain due 
to assumptions of Td and opacity compared to the 
C 18 mass, they provide independent indications 
for a lower Xco m the center of IC 342. First, 
the conversion factor estimated from dust contin- 
uum emission for GMCs B and C, are intermediate 
between the Galactic Xco value and the C 18 es- 
timates. Second, using near-infrared (NIR) spec- 
troscopy to measure the NIR extinction, Boker et 
al. (1997) obtain A„ ~ 25. From the geometry of 
the starburst (§3.1.1), it appears that GMC B is 



in front of the starburst. The NIR extinction then 
implies a column density of ~ 2 x 10 22 cm~ 2 to- 
wards the starburst region (using the Galactic con- 
version between column density and extinction). 
This value is roughly consistent with the optically 
thin C 18 value at this location. If some of the 
gas is behind the starburst then this is an underes- 
timate. But even if the starburst is in the middle 
of the cloud, the observed value is still below that 
estimated from Xco- 

Considering all the different methods, and their 
possible biases, it is concluded that XcoiJC 342) 
~ 0.7—1 x 10 20 cm~ 2 (K km s -1 ) -1 for the inner 
-75 pc of IC 342, or about 1/2-1/3 the value of 
Xco for the Galaxy. While this value is not as low 
on average as is inferred for the Galactic Center, 
it is closer to what is found towards the GMCs 
Sgr A and Sgr B2 in the Galactic Center, based 
on C 18 (Dahmen et al. 1998). 

4.3. Variation of the Xco Across the Nu- 
clear Mini-spiral 

Of equal importance to the mean Xco i n the 
center of IC 342 is its spatial variation and the 
general applicability of a constant Xco over the 
nuclear region. As noted above, there are some dif- 
ferences between the morphology of the C 18 O(l-0) 
and 12 CO(1-0) emission. C 18 in the outer arm 
regions is fainter than predicted using even the 
modified Xco suitable for the central GMCs, im- 
plying XcoiJC 342) is even lower in the diffuse 
gas. Plotted in Figure 9 (Table 5) is the estimated 
variation in XcoiJC 342) obtained from C 18 
with distance from the dynamical center for re- 
gions with detected C 18 0. Xco{IC 342) decreases 
by — 2 - 3 from the central value beyond the inner 
15 . For the diffuse gas dehind the density wave, 
C 18 O(l-0) emission is not detected at all and Xco 
may be even lower. 

Since the C 18 0(2-l) observations do not cover 
this region, we have assumed T ex = 10 K, consis- 
tent with the average seen over the mapped re- 
gion. Densities in the diffuse gas are probably 
lower than towards the GMCs, so it is possible 
that T ex for C 18 could be slightly lower than 
this due to its higher critical density. Column 
densities estimated from the optically thin C 18 
predict lower column densities for cooler temper- 
atures (even larger difference in Xco)- Moreover, 
derived N(H 2 ) are only weakly sensitive to T ex 
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for such low temperatures (Eq.(2)). Therefore, 
the possibility that C 18 is more strongly subther- 
mal in these regions cannot be the explanation for 
the differences derived in the Xco- This indicates 
that the widespread, diffuse gas seen over much 
of the mapped region is substantially brighter in 
12 CO than would be expected based on the C 18 
data. Even higher values are found in the north, 
consistent with 12 CO being especially bright in the 
north (§3.1.1). 

Dahmen ct al. (1998) find that diffuse gas in 
the Galactic Center also has a much smaller Xco 
than is estimated for Sgr A and Sgr B2, and that 
Xco is not valid across the entire region. Be- 
cause they did not have C 18 0(2-1) observations to 
go along with C 18 O(l-0) and instead relied upon 
global LVG models to predict the C 18 excitation, 
some of the structure they find probably reflects 
variations in excitation. Such variations are im- 
portant in IC 342 (§3.2.2). Hence, simply using a 
"modified" Xco is only good enough to charac- 
terize the molecular gas distribution in nuclei to a 
factor of ~2, with the amount of diffuse gas being 
overestimated. Indications of similar variations in 
the estimated Xco with distance from the center 
are also seen in M 82, and Maffei 2 (Smith et al. 
1991, Meier, Hurt & Turner, in preparation). 

4.4. An Unusually High Star Formation 
Efficiency: Dynamically Determined 
Star Formation in IC 342? 

The revised H2 column densities have implica- 
tions for the nature and efficiency of star forma- 
tion in IC 342. Since 12 CO appears to slightly 
overestimate the molecular column densities, pre- 
viously derived molecular mass fractions for the 
central region of IC 342 need to be revised down- 
ward (eg., TH92 ). Summing all the molecular 
gas (Figure 5c) over the central ~90 pc gives 
M H2 ~ 5.3 x 10 6 Mq (corrected for He and 
"missing flux"). From the nuclear rotation curve, 
the estimated dynamical mass over the same re- 
gion is M dyn ~ 5.7 x 10 7 M (TH92 ), or 
Mjjj/Mdj/n ~ 10%. The molecular mass frac- 
tion increases to M# 2 /Md yn — 30% within the 
central 45 pc. If Xco drops substantially outside 
of the central starburst region as indicated by the 
C 18 data, then the molecular arms do not make 
up a substantial percentage of the dynamical mass 
over the central arcminute. 



The dynamics (and hence the molecular gas dis- 
tribution) of the nuclear region are dominated by 
the stellar mass distribution. Recent NIR K-band 
images indicate that a small stellar bar (major 
axis ~110 pc) is present in the nucleus (Boker et 
al. 1997). The molecular gas distribution proba- 
bly reflects gas moving on barred orbits (Lo et al. 
1984; Ishizuki et al. 1990; Boker et al. 1997). 
Using the 13 CO observations, TH92 estimated 
that star formation efficiencies (SFE) towards the 
central starburst region and GMC B, in particu- 
lar, are very high. They estimate based on Ni yc 
that, at the starburst, the total mass of newly 
formed stars with masses greater than 1 M , M#, 
is ~ 10 6 M Q . With our improved H 2 mass for 
GMC B of M H2 ~ 4.0 x 10 5 M , we obtain SFE 
of M,/Mj/ 2 -2. This is an extremely high SFE. 
Alternatively, this could be an indication that the 
IMF is truncated above 1 M Q or that the starburst 
is in the process of destroying its molecular cloud. 

It is worth reemphasizing a point made by other 
studies that the sites of star formation in IC 342 
are not strictly correlated with either the global 
molecular gas column density or molecular gas 
density (Ishizuki ct al. 1990; Downes et al. 1992; 
Wright et al. 1993). Active star formation is seen 
only in two of the five central molecular clouds. 
The molecular cloud densities and temperatures 
obtained from the CO LVG analysis and HCN 
emission do not differ by more than a factor of 
^2 across the five GMCs, but the star formation 
varies substantially (§3.2.3; Downes et al. 1992). 
The locations of star formation appear to be de- 
termined predominantly by dynamics. The two 
major sites of star formation are located roughly 
perpendicular to the major axis of the small scale 
NIR bar, at the intersection with the spiral arms. 
This orientation suggests that these star forma- 
tion sites are at the intersection of the so-called 
Xi and X2 orbits (eg., Roberts, Huntley & van Al- 
bada 1979; Athanassoula 1992; Heller & Shlos- 
man 1996). 

From CO bandhead spectroscopy, Boker et al. 
(1997, 1999) find that the age of the central clus- 
ter located in the Central Trough (not to be con- 
fused with the starburst) to be 10 6 8 ~ 7 - 8 years 
old. Based on chemical enrichment models, sim- 
ilar ages are required to explain the anomalous 
C 18 abundances we find in this region (Hcnkcl 
& Mauersberger 1993). However, the sites of cur- 
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rent star formation (GMCs B & C) are probably 
younger because they are still sites of radio de- 
tected HII regions. The high emission measures 
(EM > 10 5 cm- 6 pc; Turner & Ho 1983) and the 
large implied number of O stars suggest that the 
starburst is younger than 10 Myr, perhaps signifi- 
cantly and hence younger than the central cluster. 
Star formation along the arms of the mini-spiral 
farther away from the nucleus also appears dynam- 
ically induced, perhaps by shocks associated with 
the bar orbits. The star formation along the lead- 
ing edge of the spiral arms seen in Ha emission 
can be attributed to the arm molecular gas which 
is currently upstream (TH92 ). Taking into ac- 
count the rotation speed of material in the arms, 
it is estimated that the star formation is as lit- 
tle as 1 Myr "downstream" from the molecular 
spiral arms (TH92 ). Therefore, the star forma- 
tion on the leading edges of the arms are probably 
quite young as well (and presumably roughly co- 
eval with the starburst). These result provide ad- 
ditional evidence for the following basic scenario 
for the recent star formation history of IC 342 's 
nucleus (Boker et al. 1997). Apparently a star 
formation event ^10 - 30 Myr ago formed the stars 
in the central cluster/NIR bar. The bar potential 
drove the remaining molecular gas into the nucleus 
along the x\ -orbits traced by the CO arms. The 
gas then collided with molecular gas on the per- 
pendicular x 2 -orbits, compressing the gas and trig- 
gering the younger currently seen starburst. Such 
dynamically induced star formation may then ex- 
plain why such high SFEs seen only in the central 
two intersection regions and not elsewhere. 

5. Summary 

Aperture synthesis maps of CO (1-0) and (2- 
1) in the nucleus of IC 342 are presented. This 
completes an extensive program of mapping the 
J=2-l and J=l-0 transitions of CO and its iso- 
topomers. The H 2 column density as traced by 
optically thin C 18 is similar in morphology to 
what is implied from 12 CO observations, except 
for two main differences. Firstly, the arms of the 
nuclear mini-spiral seen in C 18 (and 13 CO) emis- 
sion are more symmetric than those seen for 12 CO. 
Isotopic observations demonstrate that the asym- 
metry in brightness of the 12 CO arms are not due 
to large differences in amount of molecular gas but 
instead due to asymmetries in optically thick 12 CO 



emissivity. The northern arm is viewed with the 
starburst "lit" face towards us, appearing brighter 
in optically thick 12 CO, while for southern arm is 
viewed from the back, with the illuminated side 
away from the observer. 

Secondly, C 18 is weaker than expected from 
the 12 CO maps across the starburst region. If 
C 18 accurately traces the amount of molecular 
gas present, the standard conversion factor from 
12 CO to H 2 , Xcoi must be lower than the Galactic 
value. Comparisons between the column density 
measured from C 18 0, from 1.3 mm dust contin- 
uum emission, and from the virial theorem imply 
that X CO (IC 342) is about 1/2 - 1/3. The Galac- 
tic Xco overestimates the amount of H 2 gas in 
IC 342. For the diffuse gas component at larger 
distances from the center of IC 342, the conver- 
sion factor is about a factor of 2 - 3 lower, yet. As 
a result, a simple scaling of the Galactic conver- 
sion factor to an appropriate value for IC 342 is 
not good enough to characterize the true column 
density across the starburst. 

The mean H 2 column density for the central 
GMCs is < N H2 > ~ 4 - 5 x 10 22 cm" 2 
(A v <~ 40 — 50). This corresponds to cloud masses 
of <~ 0.5 — 1 x 10 6 M Q . Excitation temperatures 
derived from the C 18 data are T ex ~ 7 - 10 K 
over much of the central 100 pc. This value is low 
compared with what is found from the optically 
thicker 12 CO (>30 K) and 13 CO (20 K) transi- 
tions. C 18 must be subthermally excited over the 
entire range, indicating that the average densities 
of the GMCs are not higher than n c „ t (C 18 0(2- 
1)) (or ~ 10 4 cm -3 ). Densities and temperatures 
obtained for the GMCs from an LVG analysis of 
13 CO and C 18 are < n H2 > ~ 10 3 5 cm" 3 
and Tkin — 10 — 40 K. Average densities esti- 
mated from the total column densities arc consis- 
tent with these values. These values are high rel- 
ative to Galactic disk GMCs, but low compared 
to what is seen is other nearby starburst galaxies 
(eg., Jackson et al. 1995). 

Our new results indicate that star formation ef- 
ficiencies are extremely high at the starburst, with 
M*/M gas <~ 2 over the central ^45 pc. These 
measurements also indicate that molecular mass 
fractions in the center of IC 342 are lower than 
previously estimated, with M gas /Md yn ~10 % for 
the inner 90 pc rising to ~30 % over the central 
~50 pc.The locations of star formation in the cen- 
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tral ~200 pc of IC 342 appear to be determined by 
large scale dynamical forces rather than gas sur- 
face density. Star formation occurs primarily at 
the ends of the mini-spiral arms, oriented perpen- 
dicular to the large scale "bar" . This may be the 
location of the x\ — x 2 orbit crossings. Weaker 
star formation is seen downstream of the spiral 
arms apparently triggered by there passage. 

Very low 12 CO(l-0)/C 18 O(l-0) and 13 CO(l- 
0) / C 18 O(l-0) line ratios seen towards the optically 
thin Central Trough region provide evidence that 
the C 18 abundance is enhanced due to ejecta 
from the massive stars of the central starburst. 
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for reading the manuscript before submission. We 
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paper. This work was supported in part by NSF 
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Fig. 1. — C 18 O(l-0) channel map. LSR velocities 
are listed at the top right of each plane. The con- 
tours are plotted in intervals of 12 mJy bm" 1 (or 
0.12 K for the 3. 3x3. beam) corresponding to 
2a. The beam is plotted in the bottom left of the 
first plane. 

Fig. 2 — C 18 0(2-l) channel map. LSR velocities 
are listed on each plane and are the same as in 
Figure 1. The contours are plotted in intervals of 
32 mJy bm" 1 (or 0.29 K for the l."8 x l."6 beam) 
corresponding to 2a. The beam is plotted in the 
bottom left of the first plane. 

Fig. 3. — The integrated intensity of the transi- 
tions, (a) 12 CO(1-0) integrated intensity (Levine 
ct al. 1994). The 12 CO(1-0) map has been con- 
volved to the beamsizc of the C 18 O(l-0) map. 
The contours are in steps of 47.1 K km s . (b) 
C 18 O(l-0) integrated intensity. Contours arc in 
steps of 1.32 K km s _1 (2a), minus the 2a contour, 
(c) C 18 0(2-1) integrated intensity. Contours are 
the same as (b) with a 2a value of 2.68 K km s _1 . 
The large circle marks the FWHM power points 
of the mosaicked (2-1) field. The positions of the 
fitted GMCs are labeled on plot. 

Fig. 4. — The total H 2 column density maps as 
derived from C 18 (see text) overlaid on the HST 
Ha image of IC 342 (Gallagher et al, in prepara- 
tion). The contours are in steps of 6. Ox 10 21 cm~ 2 . 

Fig. 5. — Gas excitation maps. Locations with 
integrated intensities below 3a in either of the 
generating maps are blanked, (a) The C 18 0(2- 
l)/C 18 O(l-0) line ratio map. The greyscale range 
is - 1.5, with a peak value of 2.2 at the location 
of the starburst. (b) T ex derived from the C 18 0(2- 
l)/C 18 O(l-0) ratio (§3.3.1). Greyscale range is 1 
- 15 K, with a peak of 19 K. As a reference, the 
12 CO(1-0) integrated intensity map (contoured in 
steps of 125 K km s" 1 ) is overlaid on each ratio 
map. (c) The molecular gas column density de- 
rived from the C 18 O(l-0) intensity map (Figure 
3b), the T ex map in (b), and a [H 2 /C 18 0] abun- 
dance ratio of 2.94 x 10 6 . The contours are in steps 
of 2.5 x 10 20 cm~ 2 , starting with 7.5 x 10 20 cm" 2 . 

Fig. 6. — CO isotopic line ratio maps. Loca- 
tions with integrated intensities below 3<r in ei- 
ther of the generating maps are blanked, (a) 
The 12 CO(1-0)/ 13 CO(1-0) ratio map (Wright et al 



1993; MTH00). The greyscale ranges from 4 - 12 
with darker shades corresponding to higher ratios. 

(b) 12 CO(l-0)/C 18 O(l-0) line ratio map, with the 
greyscale ranging from - 75. (c) The 13 CO(l- 
0)/C 18 O(l-0) line ratio map, with a greyscale 
range of - 10. The 12 CO(1-0) contours of Fig- 
ure 6 are shown for reference. For maps that in- 
clude the 13 CO(1-0) transition, the beamsize is 
4."8 x 3. "9 (MTHOO), while the others arc those 
of C 18 O(l-0). 

Fig. 7. — Large Velocity Gradient model solutions 
for six nuclear locations. Contoured region 
represents the la confidence solution for the 
observed antenna temperature and (2-l)/(l-0) 
ratio. Solutions are displayed for both 13 CO 
(black line; MTHOO) and C 18 (grey line), with 
two assumed values of Xco/dv/dr. Abundance 
ratios of [ 13 CO/H 2 ] = 2.1 x lO" 6 and [C 18 0/H 2 ] 
= 3.4 x 10~ 7 are assumed. Solid lines are 
the solutions for dv/dr = 1 km s^ 1 pc^ 1 
(X 13CO /dv/dr = 10- 5 - 68 ; 

Xcibo/ dv/dr = 10~ 6 - 47 ), consistent with 
the data (Table 4). Dashed lines are so- 
lutions for dv/dr = 5 km s _1 pc^ 1 
(X 13CO /dv/dr = 10- 6 - 38 ; 

X C180 / dv/dr = lO" 717 ). Solutions are 
based on 4. 8 x 3. 9 resolution and a filling 
factor of 0.25 for both 13 CO and C 18 (MTHOO; 
Table 4). A decrease of a factor of two in filling 
factors will shift the solution up (<~0.4 dex) and 
to the left (<~ 10 K). Dot-dashed lines show the 
beam-averaged density estimated via the virial 
theorem, for each CMC (using the full C 18 0(2-l) 
resolution data; Table 4). 

Fig. 8. — 2.6 mm and 1.3 mm continuum maps, 
(a) 2.6 mm (110 GHz) continuum emission. The 
contours are 2 5 * (i=0,l,2...) times the 2a value of 
0.4 mJy bm^ 1 . (b) 1.3 mm (220 GHz) continuum 
emission. The contours levels are the same as 2.6 
mm, with a 2a value of 4.0 mJy bm -1 . Both maps 
are displayed of the same resolution (3 mm), (c) 
The 2 cm/ 2.6 mm spectral index map (S oc v a ). 
(d) The 2.6 mm/ 1mm spectral index map. For the 
two spectral index maps darker greyscales repre- 
sent larger values. Contours are in steps of 0.25 for 

(c) and 0.5 for (d). Negative contours are dashed. 
The maps are blanked when 2.6 mm continuum 
is less than 0.4 mJy bm -1 (c) and when 1.3 mm 
continuum is less 2.0 mJy bm -1 (d). 
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Fig. 9. — The ratio of the conversion factor 
found from C 18 observations in IC 342 versus 
the Galactic value as a function of distance from 
the dynamical center of IC 342 (Table 1). The er- 
ror bars only reflect the uncertainty in the C 18 
integrated intensity and does not include system- 
atic errors such as those associated with errors in 
derived T ex , uncertainties in the [C 18 0/H 2 ] abun- 
dance ratio or errors in the Galactic Xco- 
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Table 1 
IC 342 Basic Data 



Characteristic Value Reference 



Revised Hubble Class Sab(rs)cd 1 

Dynamical Center a(2000) = 03 h 46 m 48 s .7 ± 0. s 3 2 

(5(2000) = +68°05'46."8±2" 

2^m peak a(2000) = 03 h 46 m 48 s .3 ± 0. s 3 3 

(5(2000) = +68°05'46".8 ± 2" 

£ H ,b n 138.2°,+10.6° 1 

Vi sr 35 kms -1 2 

Adopted Distance 1.8 Mpc 4 

Inclination 30° 5 

Position Angle 39.4° 5 

B T 9.16 mag 1 



Note. — Units of right ascension are hours, minutes and seconds, and 
units of declination are degrees, arcminutes, and arcseconds. 

References. — (1) de Vaucouleurs, de Vaucouleurs & Corwin 1976; 
(2) TH92; (3) Becklin et al. 1980; (4) McCall 1989; (5) Crosthwaite, 
Turner & Ho 2000 
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Table 2 
Observational Data & 



Transition 13 


Frequency 


Tssb 


AVc/ian 




Beamsize 


Noise level 




(GHz) ' 


(K) 


(km s^ 1 ) 


(MHz) 


(arcsec; deg) 


(K / mjy BrrT 1 ) 


C 18 O(l-0) 


109.782 


400-450 


10.92 


128 


3.3 x 3.0; -43° 


0.062/6 


C 18 0(2-l) 


219.560 


430-1050 


5.46 


128 


1.8 x 1.6; -26° 


0.14/16 


3mm Cont. 


109.5 


360-430 




1000 


3.3 x 3.0; -44° 


0.004 /0.40 


1mm Cont. 


219.8 


400-1000 




1000 


1.9 x 1.6; -15° 


0.016/2 



a For observations made from 1997 October 22 - 1997 December 12 

b Phase Center #1: V L5 _R=35 km s" 1 a = 03 h 46 m 48M0 5 = +68°05'44."8 (J2000) 
#2: a = 03 ?l 46 m 49 s .01 5 = +68°05'47."8 (J2000) 
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Table 3 

Measured Intensities & Ratios 3 - 



Location 


Position 
(3 h 46"\ 68°05') 


I(C 18 (l-0)) 
(K kms- 1 ) 


I(C 18 (l-0)) 
(K kms- 1 ) 


C ls O(2-l)b 
C 18 O(l-0) 


12 CO(l-0)b 
C 18 O(l-0) 


13 CO(l-0)b 
C 18 O(l-0) 


A 


48. s 51,43."4 


8.7±1.3 


26±6.4 


1.7±0.50 


67±14 


5.8±1.2 


B 


47. s 78,46."0 


7.0±1.3 


25±6.2 


2.2±0.68 


69±16 


7.2±1.7 


C 


49. s 19,49."5 


22±3.3 


43±10 


1.2±0.34 


32±7 


5.4±1.2 


D 


48. s 93,59."2 


<4.0 


16±4 


-1.0 


>100 


>10 


E 


47. s 50,42."0 


21±3.2 


26±6.4 


0.63±0.19 


28±6 


4.1±0.90 


Arm c 


47. s 10,33."9 


4.0±1.3 






-51 


-4.7 


Off-Arm c 


48. s 10,03."8 d 


<4.0 






>110 


>7 


C. Trough c 


48. s 26,47."8 


5.3±1.3 


<8.0 




24±7 


3.2±1.0 



a The measurements of the intensities are based on the resolutions given in Table 2, while the reso- 
lutions of the ratio maps arc those of the C 18 O(l-0) (expect the 13 CO(l-0)/C 18 O(l-0) ratio which has 
a resolution of 4. 8 x 3. 9). The uncertainties reflect the larger of the absolute calibration uncertainty 
or the map noise. 

b Ratios are for integrated intensity, J T^dv, rather than the peak intensities. 
c The positions are taken from MTH00. 
d Based on a = ?> h AQ m 5 = 68°06' 



20 



Table 4 

Giant Molecular Clouds in IC 342 



Cloud 


a Q , 5 


ax b; pa 


Aui/2 




< mr N H2 > b 


Mci80° 




(03 46;68 05) 


(pc x pcf) 


(km s _1 ) 


(10 6 M Q ) 


(10 3 cm- 3 ) 


(10 6 M ) 


A 


48.51;43.4 


29 x 10; 78° 


25 


1.4 


5.8 


0.48 


Bl 


47.80;46.5 


18 x 14; 0.7° 


15 


0.46 


2.3 


0.94 


B2 


47.77;45.5 


29 x 7.8; 30° 


21 


0.86 


5.2 




CI 


49.15;51.3 


22 x 16; 17° 


27 


1.8 


5.6 


1.2 


C2 


49.19;49.5 


39 x 8.7; 180° 


26 


1.6 


5.2 




C3 


49.16;46.9 


37 x 0; 20 oa 


29 


<1.7 


<9.3 




El 


47.53;42.7 


28 x 12; 34° 


19 


0.90 


3.0 


1.0 


E2 


47.45;41.2 


19 x 16; 140° 


15 


0.54 


2.1 





a The cloud is considered unresolved since the Gaussian fit size is less than half of the beam 
minor axis of 1. 5 (13 pc). 

b The average density assumes that the cloud is an ellipsoid of line of sight dimension, £ ~ lAVab 

c Values represent the masses for the whole GMCs, since substructure is not resolved in C 18 0(1- 
0). 
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Table 5 

The Conversion Factor in IC 342 



Location 


cl8 °N H2 * 
(10 22 cm" 2 ) 


vir Na b 

(10 22 cm- 2 ) 


(lust AT C 

(10 22 cm- 2 ) 


Xco Nh 2 d 
(10 22 cm- 2 ) 


/ x MW \ 

\ Xcwo 1 


A 


3.1 


11 


-1.5 


15 


4.8 


B 


3.3 


11 


3.8 


11 


3.3 


C 


4.8 


13 f 


7.7 


16 


3.3 


D 


1.3 






13 


10 


E 


3.7 


12 




15 


4.1 


Arm 


1.1° 






7.6 


6.9 


Off Arm 


<0.55° 






3.7 


>6.7 


C. Trough 


l.l e 






5.1 


4.6 



a Using an assumed abundance of [C 16 0/C 18 0] = 250, and correcting for the 67% 
missing fiux. 

b Bascd on the sum of the masses in each subcomponent within the 12 CO beamsize 
divided by 1.36 /i# 2 (to account for He) and averaged over the beam. 

c Bascd on the dust mass, M#- 2 (Mq) estimate averaged over the beam (§3.4). 

d Using a X C o(MW) = 2 x 10 20 cm^K km s" 1 , and the 12 CO(1-0) map in Figure 3a, 
corrected for the 25% missing flux. 

°Assuming an excitation temperature of 10 K. 

f Only component C2 is included, CI and C3 are outside the 12 CO(1-0) beamsize. 
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